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Abstract 

We suggest that the spectral properties near zero virtuality of three dimensional QCD, 
follow from a Hermitean random matrix model. The exact spectral density is derived for 
this family of random matrix models both for even and odd number of fermions. New 
sum rules for the inverse powers of the eigenvalues of the Dirac operator are obtained. 
The issue of anomalies in random matrix theories is discussed. 
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In four dimensions, the spontaneous breakdown of chiral symmetry in QCD is char- 
acterized by the order parameter <qq> that relates to the spectral density of the QCD 
Dirac operator near zero virtuality, but many level spacings away from the origin [|TJ. In 
order to analyze the dynamics of the order parameter, it is therefore natural to study not 
only the asymptotic limit of the spectrum near zero virtuality, but also the way this limit 
is approached as the thermodynamic limit is taken. More concretely, we will analyze the 
spectrum near zero virtuality at finite space-time volume. 

The spectrum of the Dirac operator fluctuates over the ensemble of gauge fields. This 
raises the question whether these fluctuations may be independent of the particular dy- 
namics of the system. As we know from the study of quantum chaos and universal 
conductance fluctuations and from experimentally measured spectra of compound nuclei 
H, correlations between levels of the order of several level spacing are universal. Because 
of this they can be described by a random matrix model which has only the symmetries 
of the system as input. Our conjecture is that the same is true for the correlations of the 
eigenvalues of the Dirac operator near zero virtuality || |J [|. 

Since chirality is not defined in odd space-time dimensions, the issue of spontaneous 
breaking is more subtle in three dimensions. However, we may still consider the possibil- 
ity of a spontaneous breaking of global flavor symmetry and parity. Three dimensional 
QCD (QCD 3 ) maybe of interest to four dimensional QCD at high temperature || and 
to quantum antiferromagetism M. It may also be used to describe certain disordered 



condensed matter systems 



In this letter we would like to argue that the spontaneous breaking of flavor and/or 
parity in three dimensional QCD is related to the behavior of the spectral density near 
zero virtuality. Using the universality arguments presented above, we will derive the so 
called microscopic spectral density f|], from random matrix theory. However, in order to 
formulate the correct model we have to analyze the symmetries of the underlying theory. 
The outcome will be a Hermitean random matrix model which satisfies the 3-dimensional 
analogues of the four dimensional Leutwyler-Smilga sum rules |§. The way the anomaly 
shows up in random matrix theories will be briefly discussed. 



Consider QCD3 with Nf flavors, 



£ = - TrF 2 + J2q a (D + m a ) qi 



'a 



(1) 



a=l 



where D = ^y^D^ is the covariant derivative. The q a are two component spinors in the 
fundamental representation of SU(N C ) Color and spin indices have been suppressed. In 
the limit where all the masses are equal m a = m, there is a global U(Nf) symmetry 
of which the U(l) part relates to the baryon number, and is conserved independent of 
m a . In three dimensions, the gamma matrices are Pauli matrices. There is no distinction 
between left and right handed fermions so that U {Nf) does not have the chiral structure 
of four dimensions. 

In three dimensions, the parity operation is a reflection about one axis, say x\. This 
operation is implemented by 7072 on the spinors. As a result qq — > —qq under parity. For 
massless quarks, (ffj) is invariant under parity. However, the Dirac operator for a given 
color field is not. The total symmetry group of ([!]) is SU(N C ) x U(Nf) x Z 2 at the classical 
level. If the quark masses are arranged in pairs of opposite signs but equal magnitude, 
then (|l|) is still parity preserving. Indeed, for Nf even, a rearrangement of ([|) gives 



which is SU(N C ) x (U(N f /2) x U(N f /2)) x Z 2 at the tree level. Here, the discrete Z 2 is 
the product of the usual parity (Pi) with the interchange a —> b (E). 

At the quantum level, global anomalies can crop in. For an odd number of flavors, 
parity is either broken explicitly by quark masses, or radiatively in the massless limit, with 
the appearance of a Chern-Simons term. For an even number of flavors a rearrangement 
of the masses in doublets as in (0) causes the anomaly to vanish. In this case, parity 
is a good symmetry even at the quantum level. The parity broken phase is a screening 
phase, with free triality. The quarks are heavy and carry fractional statistics. This phase 
is related to the flux phase of quantum antiferromagnets as discussed by Wiegman ||. 
The parity symmetric phase maybe confining, with some analogy with four dimensional 



Could it be that flavor and/or parity symmetry are spontaneously broken in QCD 3 7 





QCD. 
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In so far, there is no lattice simulation to support thatQ. We suspect that for an even 
number of flavors U(Nf) x Z 2 breaks spontaneously to U(Nf/2) x U(Nf/2) x Z 2 while for 
an odd number of flavors Z 2 is radiatively broken by the anomaly in the massless case, 



and explicitly in the massive case. By a generalization of the Vafa-Witten theorem [13 



to QCD 3 it follows that the absolute values of all condensates are equal. As follows from 
the derivation of the Banks-Casher relation, the quark mass and the condensate for each 
flavor have the same sign if the symmetry is broken spontaneously. Therefore, we can 
introduce the order parameter 

I N f 

\n = w J2\<^q a >\. (3) 

/ a=l 

Throughout, we will study the spectrum of the Dirac operator at finite space-time 
volume. In the chiral limit the spectrum is well defined and contains information on the 
structure of the vacuum, including the way the flavor symmetry is broken. In the large 
N c limit, a rerun of the Coleman- Witten argument |14[] for QCD% suggests that in the 
even case the symmetry is broken according to U(Nf) — > U(Nf/2) x U(Nf/2). This is 
confirmed below using an effective Lagrangian method in the saddle point approximation. 
In this spirit and for even Nf, the condensates will be combined in pairs of opposite sign. 

To understand the generic behaviour of the Dirac spectrum near zero virtuality, we 
note that the massless Euclidean Dirac operator iZ)[A] in an arbitrary background field 
is Hermitean and does not commute with PiE. By analogy with four dimensions, the 
random matrix model with the symmetries of (|2]) is 

Z(m) = f VTP(T) [] det(*T + m a ), (4) 

J a=l 

where the Haar measure is over N x N hermitean matrices, and the mj are the eigenvalues 
of the mass matrix. The weight distribution P(T) is chosen to be Gaussian consistent 
with no additional input but the symmetries of the system. We note that the fermion 
determinants make the integrand not necessarily positive. This ensemble with the deter- 
minant replaced by its absolute value is also known as the generalized gaussian ensemble 

m 



1 However, lattice calculation have been performed for QED$ H], and a nonzero value of (qq) was 
found. 
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The order parameter (|j) follows from 



1 d 

E = — lim lim — log-Z" (5) 

m a ^0N^oo N dm a 

If we were to define the continuum "spectral density" by 

p c (A)= hm hmd$>(A-A n )> (6) 

m— >0AT— >oo N 

1 ' n 

where the expectation value is over the partition function (|2]), then (|5|) can be rewritten 

as 

E = iirpciO) - P / dX^- (7) 

For even spectra, the principle value part vanishes. 

The partition function can be evaluated if we recall that the general decomposition 
for Hermitian matrices is T = UAU^ , where A is a diagonal matrix. Using the eigenvalues 
and eigenangles of T as new integration variables, the partition function can be rewritten 
as 

r Nf NT 2 
Z(m) = Y[ dX kE\h-Xi\ 2 UIK iX k + m a )ex P ( — £ X 2 k ) (8) 

k k<l k a=l ^ k 

The one-point function p(X) is defined as the integral over all eigenvalues in this partition 
function except one. Its normalization can be expressed through Z = J dXp(X). From 
(§) it follows immediately that p{—X) = (— 1) fp(X). For an even number of flavors 
p(X) is positive in the chiral limit and can be interpreted as the spectral density of the 
Dirac operator. For an odd number of flavors p(X) is still an even function for even N 
but is not necessarily positive definite. For odd iV and an odd number of flavors it is 
an odd function. In this case Z(0) = and Z{m) ~ d m Z\ m=Q . In the calculation of 
the condensate (see ([5[)) the derivative of the partition function cancels and to leading 
order in m we obtain E = lim TO ^o limjv->oo( — 1)/Nm = 0. Thus, chiral symmetry remains 
unbroken in this limit. 

Before discussing the spectral density related to (|]), we will derive the finite volume 
partition function in the static limit, as for QCD4 0. This is achieved by rewriting the 
fermion determinant as Grassmann integrals, and averaging over T. The result is a four- 
fermion interaction. The effective partition function is obtained after bosonisation and 



a saddle point approximation. To leading order in 1/JV only the saddle point with the 



eigenvalues of a in opposite pairs contributes |L5[ . The resulting partition function has a 
'hyperbolic' symmetry and is given by 

Z(M)= [ DUexp(NXTrMUIlP), (9) 

JU£SU(N f ) 

where we have extended the integral over the coset U(Nf)/U(Nf/2)xU(Nf/2) to SU(Nf), 
and I = diag(ljVf/2; — ljVy/2)- Equivalently, the result (||) could also be arrived at using 
general symmetry arguments ||. 

The linear term in M. in the expansion of Z(M) vanishes. The integrals over SU(Nf) 
that occur in the terms of 0(M. 2 ) are well known. We find 



Z(M) = Z{0) 



Y?N 2 N 2 (I k „ 1 „ \ 

1 H — — TiM Tr 2 .M + 

+ 2 N]-l \N f Nf J + 

On the other hand, the QCD% partition function can be expanded as 



(10) 



Z(m) = Z(0)( 1 + *Tr.M £ 1 + i(Tr.M 2 - TrlM) £ i - ^Tr 2 ^ £ J- + •••]), (11) 
resulting in the sum rules 



Note that the average over the second term in ([0]) vanishes, in agreement with the 
effective partition function. 

For odd Nf it is not possible to organize the saddle points in opposite pairs. Because of 
the Jacobian and to leading order in 1/N the saddle points occur at (Nf + l)/2 eigenvalues 
with ±1/S and (Nf — l)/2 eigenvalues with =Fl/S. Both sets of saddle points cannot 
be transformed into each other by a unitary transformation. However, all saddle point 
in one set are connected by a unitary transformation. The sum of the two saddle points 
leads to the partition function 

Z(M)= [ DU cosh(NT,TrM Ulrf),, (13) 

JueSU(N f ) 

for even N, and cosh — > sinh for odd N. The matrix I is now diag(l(jv>+i)/2 3 — l(A r / -i)/2)- 
The finite volume partition function is therefore based on the coset U(Nf) /U((Nf — l)/2) x 
U((Nf + l)/2). Sum rules follow from comparison to the QCD% partition function 

^?AT^ md ^§AA = "^ <14) 



for even values of N. 

With these results in mind, we now turn to the evaluation of the spectral density using 
random matrix theory. The microscopic spectral density and its correlation functions will 
be the master formulae for all sum rules. The construction of the spectral density, can be 
achieved with the help of the orthogonal polynomial method from random matrix theory 
(see for example 0). For even Nf we proceed as for the Gaussian Unitary Ensemble 
The result for the spectral density can be written as 

N-l -I 

p(X) = E -Pk(X)Pk(X)X Nf exp(-a 2 A 2 ), (15) 
fc=o rfc 

where the P k axe orthogonal polynomials that satisfy 

dXP k (X)Pi(X)X N f exp(-a 2 A 2 ) = r k 5 kh (16) 



and r k is a normalization factor. For even Nf these polynomials can be expressed in terms 
of the generalized Laguerre polynomials 



Nf-l JV / + 1 



P 2k (X) = L k 2 (a 2 A 2 ) and P 2fc+1 (A) = XL k 2 (a 2 X 2 ). (17) 
Using (|T6| ) and (p!7|) the normalization factors are determined readily. 



The spectral density follows immediately from (15). It can be written as the sum 
of two terms. Both sums can be evaluated with the help of the Christoffel-Darboux 
formula. The result is straightforward, and yields an exact analytical expression for the 
level density. Its microscopic limit 

*w=&^- (18 > 

follows from the asymptotic properties of the generalized Laguerre polynomials. The 
result for even Nf is |12 



Ps{z) 



Z 2 z 



J% f -l(z) — J nj+i (z)J n s -z (z) + J% f +i(z) — j N f +3 (z) J N f -i (z) 



(19) 



It reproduces the first sum rule of ([12]). 

For Nf odd it is not possible to construct a set of orthogonal polynomials satisfying 
(|16|). For example, the zeroth order polynomial cannot be orthogonal to the first order 
polynomial. Therefore, in the Vandermonde determinant, we replace the 2k — l'th row 



by the difference of the 2fc'th and the 2k — l'th row and the 2/c'th row by the sum of 
the 2k — l'th and the 2fc'th row. The polynomials in the 2k — l'th and the 2/c'th row 
therefore have the same degree. In this case the level density is again given by (|15|) but 
the orthogonal polynomials P 2 k-i and P 2k have the same degree. The construction of 
these polynomials is straightforward 

Nf Nf 

Pik-i = (1 - Aa)L & 2 (AV) and P 2k = (1 + Xa)L k 2 (AV). (20) 
Repeating the steps leading (|19"D we arrive at the microscopic spectral density 

p s (z) = (-1) N/2 ^- (j 2 Nf/2 (z) - J Nf/2+1 (z)J Nf/2 -i(z)) , (21) 



which reproduces the first spectral sum rule of ([L4]). The case of odd iV and odd Nf will 
be discussed elsewhere. 

Finally, we would like to comment on the issue of the anomaly. In three dimensions, 
det(zZ?3[v4]) is noninvariant under an adiabatic switch-on or -off of a large gauge trans- 
formation |16|]. The gauge noninvariance is usually followed by a spectral flow in the 
fermionic spectrum with level crossing at zero. This means an overall change in the sign 
of the fermion determinant. Thus the anomaly and the phase of the fermion determinant 
are related. Indeed, consider the case when one eigenvalue say A, crosses 0. The relevant 
factor in our case that determines the change in phase is given by ]T^2i(A + im a ). When 
the crossing is completed with \m a \ <C |A|, the determinant acquires the extra phase 

Nf 

exp(±z7r sgnm a ). (22) 

a=l 

This result is to be contrasted with the continuum result [ |T6| 

N f 

exp(±zvr^sgnm a W CS [A]), (23) 

a=l 

where W CS [A] is the Chern-Simons action in Euclidean space. The imaginary character 



of ( |23"D in Euclidean space, follows from the fact that W cs is T-odd. The sign ambiguity 
in our case, corresponds to the sign ambiguity left by the Pauli-Villars regulator in the 
continuum. Thus, the net sign effect in random matrix theory amounts to a Chern-Simons 
term. 
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To conclude: we have argued that on the basis of universality arguments that the 
spectral density of QCD 3 near zero virtuality follows from a hermitean random matrix 
model. We have explicitly constructed the microscopic spectral densities for an even and 
odd number of flavors. The resulting sum rules are in agreement with the expected results 
following from an effective Lagrangian formulation based solely on symmetry arguments. 
For an even number of flavors, we have suggested using the effective Lagrangian formula- 
tion, that U(Nf) is likely to be maximally broken to U(Nf/2) x U(Nf/2). We have also 
shown that the Chern-Simons term has a natural explanation in the context of random 
matrix theory. 

Finally, we want to note that the random matrix model used in this work is based on 
the assumption that flavor symmetry is broken spontaneously. It is this assumption that 
leads to the finite volume static partition function as quoted in the text. Since, in general, 
odd dimensions do not sustain semiclassical physics, this leads us to the questions: what 
is the mechanism behind the spontaneous breaking of flavor symmetry? Could it be that 
this mechanism is also responsible for the spontaneous breaking of chiral symmetry in 
four dimensions ? These points deserve further investigation. 
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